Josephson array of mesoscopic objects. Modulation of system properties through the 

chemical potential 
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The phase diagram of a two-dimensional Josephson array of mesoscopic objects (superconducting 
(3J[)' granules, superfluid helium in a porous medium, traps with Bose-condensed atoms, etc.) is examined. 

Quantum fluctuations in both the modulus and phase of the superconducting order parameter are 
taken into account within a lattice boson Hubbard model. Modulating the average occupation 
number no of the sites in the system (the "number of Cooper pairs" per granule, the number of 
atoms in a trap, etc.) leads to changes in the state of the array, and the character of these changes 
depends significantly on the region of the phase diagram being examined. In the region where there 
are large quantum fluctuations in the phase of the superconducting order parameter, variation of 
the chemical potential causes oscillations with alternating superconducting (superfluid) and normal 
states of the array. On the other hand, in the region where the bosons interact weakly, the properties 
, of the system depend monotonically on no. Lowering the temperature and increasing the particle 

interaction force lead to a reduction in the width of the region of variation in no within which the 
system properties depend weakly on the average occupation number. The phase diagram of the 
array is obtained by mapping this quantum system onto a classical two-dimensional XY model 
with a renormalized Josephson coupling constant and is consistent with our quantum Path-Integral 
Monte Carlo calculations. 

B 

^ i 1. INTRODUCTION 

Advances in micro-lithographic techniques have made it possible to create regular arrays of mesoscopic Josephson 
objects. These systems, which are under active experimental and theoretical study, include superfluid helium in 
porous media 1 ' (see Ref. 2 and the literature cited there), lattices of mesoscopic Josephson contacts, 3 ' 4 and ultrasmall 
superconducting granules. 5,6 One interesting physical realization of a Josephson array consists of Josephson junctions 
in a structure with superfluid 3 He created by means of lithography. 7 Major advances in experiments with Bose 
condensates of atoms cooled by laser irradiation followed by evaporation 8-10 suggest that it may be possible to 
^—1 j fabricate a Josephson array of close magneto-optical traps with Bose-condensed atoms 2 - 1 or with clusters of Bose- 
condensed atoms cooled and localized at the nodes of a system of standing electromagnetic waves. Finally, another 
remarkable realization of a Josephson array might be a system of Josephson coupled "lakes" of Bose-condensed excitons 
in single or double quantum wells located at the minima of the random field created by the roughness of the well 
surfaces, i.e., in "natural" quantum dots, 12 or in an array of artificial quantum dots. 

For concreteness, the following discussion considers the example of a system of superconducting mesoscopic granules 
or Josephson junctions, but the results apply to all the systems mentioned above. We consider a regular array of 
mesoscopic granules situated on a conducting substrate and separated from it by a thin dielectric layer. A voltage 
applied to the conducting substrate serves as the chemical potential of the Cooper pairs, which determines the 
average occupation number no of the granules in the system. 13,14 For example, in the case of superfluid helium in 
a porous medium, the chemical potential of the atoms can be varied (because of the contribution of the van der 
J-^ ' Waals interaction) by changing the thickness of the layer of adsorbed helium. 15,16 The state of an array of mesoscopic 
^ , traps with cooled atoms can be controlled by changing the average number of atoms in the system (by, for example, 
capturing them from an external flow). 

The systems being examined here, in which the character of the processes taking place is determined by the boson 
degrees of freedom, are conveniently described by a lattice boson Hubbard model 17,18 with the Hamiltonian 

£ JJ 2 

H = - (2a^i(ii — aUdj - a^-a,-) + — ^ (a^a,) — M ^ a ^i a i (!) 

<i,j> i i 

In this model site i corresponds to one superconducting granule or pore with helium, to a single trap with a Bose 
condensate, etc. The operators a\ (of) are the Bose creation (annihilation) operators of an "effective" boson at site 
z = l, ./V 2 of an A x A lattice. The first term in the Hamiltonian takes into account the kinetic energy of the particles, 
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which corresponds to the Joscphson tunneling energy and is described by the parameter t. The second term in Eq. 
(1) describes the interaction of effective bosons at a granule with a characteristic energy Ui > 0. 

The model (1) is interesting in that it can be used to study the properties of arrays of mesoscopic structures in 
which the relative fluctuations in the modulus of the superconducting order parameter are large. In this regard, we 
note that the quantum XY model is justified only when these fluctuations are small, 14 i.e., in the case of arrays of 
macroscopic granules. 

A lattice system (1) with Mott- insulator and superconducting phases at T = has been studied both 
analytically 17-19 and by computer simulation. 20,21 In this paper we shall be interested in the properties of the system 
(1) at finite temperatures, digressing from the interesting phase transitions at T = 0. 22 In our earlier papers, 23 the 
investigation was limited to the case of integer (commensurable) populating, in which the average number of bosons 
per site (granule), n = {a\a,i), is a whole number. At T = 0, adding even a single particle to an arbitrarily large 
system changes its properties in a fundamental way. Specifically, a system with a noninteger average number of 
bosons per granule remains superconducting for arbitrary values of U jt for the interaction between the particles. 17-22 
It is evident that this surprising behavior should be present even in the limit of T = 0. In fact, as will be shown 
below, at finite temperatures the properties of the system vary little within an interval no = k ± 8no about an integer 
population, whose width 2Sng decreases as the temperature is lowered. 

The purpose of this paper is to study the changes in the character of the ordering in an array of granules as 
the substrate voltage (chemical potential of the effective bosons) is varied. Here we shall not use the simplifying 
assumption of small relative fluctuations in the modulus of the superconducting order parameter. The results given 
below correspond to an array of mesoscopic objects for which the root-mean-square fluctuations in the number of 
particles are comparable to their average number. In Sec. 2 we present results from mean-field calculations. The 
method used there corresponds to mapping the initial boson model (1) onto an effective classical XY model with a 
renormalizcd Joscphson coupling constant. To refine the results of the analytical calculations, we use the quantum 
Path-Integral Monte Carlo method (see Sec. 3). A discussion and comparison of the results in Sec. 4 completes the 
presentation. 



2. BOSON HUBBARD MODEL IN THE MEAN-FIELD APPROXIMATION 



A qualitative approximation for the phase diagram of the model (1) can be obtained using an approach described 
previously. 19,23 ' 24 In terms of this model, the boundary of the ordered state is located where the local density of 
the superconducting component vanishes in the effective functional describing long-wavelength excitations of the 
system. The latter can be obtained in the usual way using the Hubbard-Stratonovich transform 25 followed by an 
expansion of the effective functional in components of the fluctuating field. 19,24 The condition that the local density 
of the superconducting component vanish for the system described by the Hamiltonian (1) yields an equation for 
determining the boundary of the ordered state: 



4^ exp|-0.5<z 2 (n-r)) 2 /f I - cxp j -0.5<f(n + 1 - ?)) 2 /T \ {n + l)/(2(n - fj) + 1) 



f] nexp\-0.5q 2 (n-fj) 2 /f) 



(2) 



n=0 



Here we have used the following independent dimcnsionlcss parameters, which specify the state of the system: 



k b T 

t ' 



a 2t 



q = y/U/t, T 

The average number of particles in the granules of the system in the disordered state is given by the equation 

J2 nexp \ -0.5<7 2 (n - fj) 2 /T \ 



n = 



£ exp{-0.5g 2 (n-ry) 2 /f| 

n=0 > 



(3) 



The solid curves in Fig. la represent the phase diagram of the system (1) in the variables {^JU/t, u/U} obtained by 
solving Eq. (2) for different temperatures ksT/t. The region of low values of U/t (low particle interaction energies) 
corresponds to the superconducting (S) state of the system. 
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FIG. 1. Phase diagram for the Hubbard model (1) in the coordinates {yJU/t, n/U}. S — the superconducting state; 
N — the normal (metallic) state; I — Mott insulator (hatched region in Fig. b). (a) Mean- field calculation. The solid 
curves were obtained by solving Eq. (2) and correspond to vanishing of the local density of the superfluid component. The 
Kosterlitz-Thouless topological transition (6) takes place on the dashed curves, (b) Monte Carlo calculations for feT/t = 0.8. 
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FIG. 2. Phase diagram for the Hubbard model (1) in the coordinates {q, n/U} and {q, Q}. (a) Mean-field theory calculation. 
The solid curves were obtained by solving Eq. (2). The dotted curves correspond to the joint solution of the system of Eqs. 
(2) and (3). (b) Monte Carlo calculations. The filled squares are the phase diagram for the 2 + 1-dimensional XY model (for 
integer values of n = k S> l). 34 The system with {U /t, ksT/t} = {3.5, 0.8} and a variable chemical potential fi/U moves along 
the dotted curve (see Figs. 4 and 5 below). 



As k B T/t — > 0, the disordered state of the system corresponds to an integer average number of particles per 
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granule, no = k, and determines the domain of existence of the Mott insulator (I). 17 In this limit, the superfluid 
state of the system corresponds to the case of incommensurate populations, i.e., to a noninteger average number of 
particles per granule. Figure la shows that for T = and half-integer values of the chemical potential fi = 0.5 + k, a 
superconducting state exists for arbitrarily strong interparticle interactions, in accord with earlier works. 17,18 

At finite temperatures an increase in the boson interaction force leads the system into a disordered state for arbitrary 
values of the chemical potential (see Fig. la). It is also clear from the figure that as the temperature is increased, the 
domain of existence of the ordered state is shifted toward higher values of the chemical potential. 

The transition to the case of a system of macroscopic granules corresponds to increasing the particle density no and 
reducing the role of the fluctuations in the modulus of the order parameter. It is most convenient to follow the changes 
in system (1) with increasing no in the {q, 0} plane, where we use the dimcnsionless temperature 8 = ksT jtno and 
the quantum parameter q = \JUjtrio, which are the control parameters that also determine the state of the quantum 
XY model. The corresponding phase diagram is shown in Fig. 2a. As can be seen from this figure, for any values of 
U the estimate of the boundary of the ordered state in the Hubbard model according to mean-field theory lies above 
the corresponding limit in the XY model and approaches it as the average population no of the lattice sites increases. 
Our calculations confirm that the phase diagram is periodic in the parameter fi/U when no 3> l. 14 

The dotted lines in Fig. 2a comprise a family of curves, at whose points a system with density no becomes disordered. 
The points where they intersect the O = const lines determine the phase diagram of a system in the {q, 0} plane 
corresponding to n particles per granule. A similar analysis shows that for an incommensurate boson density at 
low temperatures (see below), an ordered (superconducting) state of an array exists for arbitrarily large values of q, 
i.e., for arbitrarily large quantum fluctuations in the phase of the superconducting order parameter in terms of the 
quantum XY model. 

This approach only yields a qualitative estimate of the characteristic features of the phase diagram of this system. 
A comparison with the results of a numerical simulation (see below) shows that Eqs. (2) and (3) give a greatly 
overestimated value for the disorder temperature C (<7; v/U). In order to obtain more accurate quantitative estimates, 
it is necessary to determine the temperature at which the global (rather than local, as in the method described above) 
superfluid density of the array vanishes. This temperature can be estimated as the Kosterlitz-Thouless temperature 
for the topological phase transition according to the classical XY model, onto which the initial system is mapped by 
expanding the effective Ginzburg-Landau functional for weak fluctuations in the phase of the order parameter. 

Using the approach in Refs. 19, 24, and 26 it is easy to show that the effective action of the classical two-dimensional 
XY model sought is given by 

/ tA 2 
S(W) = J 5 L 'E\k\ 2 <Pj;<P-i*>JxY E (l-coB(<p i -<p j )),J XY WU;t/U;k b T/t) = —, (4) 

k <i,j> 

where Jxy is the coupling constant for the effective XY model, which depends on the dimensionlcss parameters n/U, 
t/U, and ksT/t. The local superfluid density A 2 /4 of the system is described by the relation 

tr { (fit + o)e-^-/} u 

A = — ^ : — , L ,H m} = -n 2 + (2t - /i)h - iA(fit + fi) (5) 

tr I e -P H ™t \ 2 

A plot of the Kosterlitz-Thouless phase-transition temperatures in the effective XY model 27 specifies the boundary 
sought of the superfluid state for this array of granules: 

k b T c = Q.mJxY^/U- t/U; k b T c /t) (6) 

Estimates given by Eqs. (4) and (5) are shown as dashed curves in Fig. la. Note that, although the two approaches 
discussed above yield a similar qualitative behavior of the boundary of the ordered phase, the temperature of the 
topological phase transition in the effective XY model (6) is considerably lower than the temperature at which the 
local superfluid density vanishes according to Eq. (2). A comparison of the phase diagram obtained in this way with 
the results of Monte Carlo calculations (see below and Fig. 1) shows that they are in fair quantitative agreement. 



3. QUANTUM MONTE CARLO METHOD. MEASURABLE QUANTITIES 

The Trotter discretization procedure makes it possible to estimate all the thermodynamic averages of the observables 
of a D-dimensional quantum system in terms of a classical D + 1-dimensional system, where the product of the matrix 
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elements (calculated approximately) of the high temperature density matrix serves as the Boltzmann weight of the 
configurations of the corresponding classical system. For studying the properties of the model (1), we shall use the 
"checkerboard version" of the quantum Monte Carlo method. (A detailed discussion of the discretization procedure 
and the organization of the Monte Carlo step during the simulation of systems of lattice bosons in a large canonical 
ensemble is given elsewhere. 28 ) In this method, the degrees of freedom of the discretized system are the occupation 
numbers {nf } of the sites of the N x N x 4P three-dimensional lattice formed by 4P-fold multiplication of the initial 
N x N lattice along the imaginary time axis. The number of subdivisions P was chosen so that the parameter 
e = q 2 /P 2 <d 2 , which characterizes the discretization error, would be less than 0.06. 

The density v s of the supcrfluid component was calculated at each computational point of the phase diagram, whose 
position is specified by the parameters {y/U ft, ksT/t} and the chemical potential \xjTJ . To find this quantity, we 
used both the fluctuations in the topological winding number 20,28 and the correlation function of the paramagnetic 
current. 29 We found that when the average particle density at the boundary no < 2 and q > 2, the statistical errors in 
the second method were considerably higher than the errors in determining the supcrfluid density from the fluctuations 
in the winding number, rendering it unsuitable. 

We also measured no, which is controlled by the chemical potential of the system, and the compressibility modulus 
k, which is defined as 

n = hTdno/d^ = — (EE ) - M 2 ( ? ) 

\ p=o i I 

It turned out to be convenient to make the measurements with U/t = const and fcsT/t = const and variation of the 
chemical potential fi/U. As an example, in Fig. 2b the line along which the system moves for {U/t, ksT/t} = {3.5, 0.8} 
is plotted in the coordinates {q, 0}. The location of the system on this line for a given value of the chemical potential 
can be determined by measuring the average number no of particles per granule. In addition, we performed a number 
of calculations for fixed no (in a canonical ensemble). It might be expected that for the same particle density no the 
results of the simulation would be independent of the choice of the ensemble for sufficiently large systems. We tested 
this assumption and found for a system of dimension N x N = 6 x 6 that the difference in the measured quantities is 
less than 10% in the region of variation of the control parameters of interest to us. Therefore, in analyzing the results, 
data obtained by assigning the density n (within a canonical distribution) and by assigning the chemical potential 
(which corresponds to using a grand canonical ensemble) can be used simultaneously. 




FIG. 3. Plots of the fraction v s of the superfluid component as a function of the average occupation number no. The 
solid curves were obtained by interpolation of the data with a fourth-order polynomial. The lengths of the horizontal arrows 
correspond to estimates of Sn(q, O) = n{q, 0)|„ o= i. See Eq. (8). 
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4. DESCRIPTION AND DISCUSSION OF RESULTS 



We first consider the computational results for q < 1.5, where the interparticle interaction is weak and mean-field 
theory (see Figs, la and 2a) predicts a monotonic dependence of the phase-transition temperature on the particle 
density no. Figure 3 shows the calculated density of the superfluid component as a function of average occupation 
number, v s (n ), for q = 1 and 9 = 1 (unfilled symbols). As the occupation number increases, v s approaches a constant 
equal to the helicity modulus -f(q, 9) in the 2 + 1-dimcnsional (quantum) XY model. 30,31 We found previously 23 that 
in this region of the phase diagram, this limit, which corresponds to a small contribution of the quantum fluctuations 
to the modulus of the order parameter, is approached already when n = 4-5. The monotonicity of iAs(n ) suggests 
that the results of experiments on a system of isolated granules will not differ greatly from those on a system of 
granules on a substrate with an applied potential. It turned out that the system behaves similarly up to q ~ 2.3. 
As the quantum parameter q is raised further and the temperature 9 is lowered, v s (no) ceases to be a monotonic 
function of the average occupation number uq. Characteristic oscillations in v s {no) with minima at integer values of 
n = k are noticeable in Fig. 3 ({q, 9} = {2.75,0.5}, filled symbols). For sufficiently high boson densities (n > 7 
in the region {q, 9} = {2.5, 0.5}; see Ref. 23) there is a transition to the quasiclassical limit and the density of the 
superconducting component v s (no) becomes a periodic function of the average occupation number with a period of 
unity. 14 

Figure 4a shows plots of the superfluid density as a function of the chemical potential of the system, which is 
proportional to the voltage applied to the substrate. The squares refer to a system with \fUjt = 2, ksT/t = 0.8 and 
the triangles to \fU~Ji = 2.5, ksT/t = 0.8. The feature at n ~ 1 [a dip on the plots of v s {^/U) and k(/i/U) for at 
n/U ~ 1.2] is poorly seen for q ~ 2 (squares in Fig. 4) but becomes clearly evident for q ~ 2.5 (triangles). The figure 
confirms that the case of commensurate populations corresponds to lower densities of the superfluid component, i.e., 
increasing the deviation of the average boson density from integer values of n = k leads to spreading of the phase 
diagrams for the model (1) in the {q, 9} plane. 

As the quantum parameter q is increased further, the differences in the properties of the system for integer and 
noninteger boson densities become increasingly more significant. In fact, based on the results from mean-field calcu- 
lations (see Fig. la), we can assume that there is a value of U/t for which the \fU~ft — const line intersects the region 
of the disordered state with no — 1. Further increases in the interaction constant U should lead to the possibility of 
an intersection with the disordered region having n ~ 2, etc. The computational results shown in Fig. 5 confirm 
this assumption. Note that for \/TTJt = 5.0 (triangles in Fig. 5) and integer values of n/U, changes in the chemical 
potential lead to essentially no change in the average number of particles at the sites in the system. This feature, 
which is characteristic of an insulator, can also be seen in Fig. 5b, which shows a plot of the compressibility modulus 
Kj as a function of the chemical potential (J,/U. 

At T — a lattice boson system without disorder undergoes a superconductor- insulator phase transition. 17 As 
for the system being studied here, one can assume that as the temperature is raised, the domain of existence of the 
Mott insulator decreases and shifts toward larger U/t, so that, for example, on moving along the [i/U = k line the 
superfluid phase is replaced by the normal (metallic) phase. Further increases in the interaction constant, under 
which the magnitude of the Mott dielectric gap increases so much that thermal excitations become unimportant, lead 
to crossover formation of the insulator state. 32 An estimate of the location of the boundary of the insulator state 
for no = 1, based on our calculations, is shown in Fig. lb. Several constant-density (no = const) contours are also 
shown for comparison. A more detailed study of crossover-type metal-insulator transitions requires that systems of 
substantially larger size be examined. 24,33,34 

An analysis of Figs. 1-5 shows that interesting effects caused by incommensurate populating of the sites in the 
system occur only in quite strongly interacting systems and at sufficiently low temperatures. For the boson Hubbard 
model (1), the domain in which they exist can be estimated by the inequalities q > 2.3 and 9 < 0.7. Quantitatively, 
the magnitude of the deviation Sn of the average particle density from an integer value of n — k at which a significant 
change in the system properties will be observed can be expressed in terms of the compressibility modulus: 

hT = «|„ o=fc (8) 

no — k 

It is known that the compressibility modulus [which is inversely proportional to the phase fluctuations in the 2 + 1- 
dimcnsional XY model; see Eq. (7)] falls off substantially as the quantum parameter q becomes greater and the 
temperature 9 is reduced. We found that the following estimates hold for the Hubbard model (1) when n = 1: 
Sn w 1.2 for {g,9} w {0.5,1}, Sn « 0.3 for {q,G} « {2.5,0.8}, and Sn w 0.06 for {q, 9} « {2.75,0.5}. As 
an illustration, Fig. 3 shows the values of Sn(q,Q) found for the points {q, 9} = {1,1} (unfilled symbols) and 
{q, 9} = {2.5, 0.8} (filled symbols). The figure demonstrates the fair agreement between the theoretical estimate (8) 
and the Monte Carlo calculations. 
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FIG. 4. (a) The superfluid density v s as a function of the chemical potential fi/U. (b) Average particle density no 
(filled symbols) and compressibility modulus k (unfilled symbols) as functions of the chemical potential n/U: squares - 
{U/t,ksT/t} — {2.0,0.8}, triangles — {U/t,kBT/t} = {2.5,0.8}. Spline interpolations are shown for visual convenience. 
When not indicated the statistical errors are smaller than the sizes of the corresponding symbols. 
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FIG. 5. (a) The superfluid density v s as a function of the chemical potential fi/U. (b) Average particle density no 
(filled symbols) and compressibility modulus k (unfilled symbols) as functions of the chemical potential (J./U: squares - 
{U/t,k B T/t} = {3.5,0.8}, triangles — {U/t, k B T/t} = {5.0,0.8}. 

The above results yield the phase diagram of the system shown in Fig. lb (in the coordinates {\/U/t, /J./T} for 
ksT/t = 0.8) and in Fig. 2b (in the coordinates {q 7 6}). The location of the boundary of the ordered superconducting 
state was estimated from the universal jump in the superfluid density 30 and from the location of the peak in its 
temperature derivative. 35 



9 




0.0 0.5 1.0 1.5 



FIG. 6. Superfluid density of the system as a function of O = fcsT/tno for q = 3.0; unfilled symbols — no = 1.306; filled 
symbols — n = 1.194. The dotted line indicates the straight line v s = 2Q/n. The inset shows the calculated superfluid density 
of the system as a function of fe^T/i for WU/t = 2.5 and n/U = 0.75. 

Quantum phase transitions (at T = 0) in the two-dimensional Hubbard model (1) are determined by the critical 
properties of the corresponding effective three-dimensional system. 17 ' 32 At finite temperatures a Josephson array will 
display Kostcrlitz-Thouless critical behavior on some c (/i/£7, q) curve. We now try to evaluate the effect of quantum 
fluctuations on the temperature of this transition, assuming it to be rather low (see below). There is interest in two 
cases: (a) the system has a commensurate particle density n — k and at some temperature @ c {q) it undergoes a 
transition from the superconducting to the normal state, and (b) a phase transition takes place at 6 c (no) because of 
a change in the density no. 

For the region near the point q* Y « 2.5 of the quantum transition in the 2 + 1-dimensional XY model (see 
Fig. 2b), the temperature C (</) of the topological Kosterlitz-Thouless phase transition has been estimated 24 as 
©c(<z) ~ \q — Q^l'' 1 W1 th Ci ~ 0.67. This estimate was derived under the assumption that the system temperature 
is less than the temperature of the 2D — > 3D crossover; i.e., O < &sd, where ©3D ~ \q — q* Y \ u with v w Ci- This 
makes it possible to estimate the density of the superfluid component, v s (q, 0), which determines the phase-transition 
temperature c (g), as v a (q, 0) w v s (q,0). Evidently, similar arguments apply in the neighborhood of the points 
q^\n =k of the quantum phase transitions in the Hubbard model (1) for integer populating of the array granules, i.e., 
n = k (q^\n =i ~ 2.8; see Ref. 20). Thus, we may expect that 

^k)-\ q c H \ no=k - q \°- 67 (9) 

Similar arguments can also be applied to the case of incommensurate boson densities, n ^ k. It has been shown 17,20 
that v s ~ \n — k\& with ( 2 ~ 1-0 for q > q^ . Thus, the following relation holds for the temperature c (no) of the 
topological Kostcrlitz-Thouless phase transition: 

c (n o )~|n-fc| 1 - (10) 

Our quantum calculations are in qualitative agreement with the predictions of Eqs. (9) and (10), but it is difficult to 
confirm their validity with sufficient accuracy because of the large errors in determining the position of the quantum 
phase-transition line fj,(U/t;T = 0). 

There is great interest in the question of the existence of reentrant superconductivity, for which, within some range 
of variation of the quantum parameter q, disorder sets in not only as the temperature is raised, but also as it is 
lowered. The existence of reentrant effects has been predicted a number of times within the quantum XY model (see 
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Ref. 13 and the literature cited therein), but as far as we know, computer simulation cannot unequivocally confirm 31 ' 36 
or refute 3 the existence of this phenomenon. The earlier numerical calculations of the Hubbard model did not reveal 
low-temperature instability or reentrant superconductivity. 20 ' 23 For this paper we studied the low temperature region 
q w 2.5, 6 < 0.5 of the boson Hubbard model both for different noninteger occupation numbers and for fixed values 
of the chemical potential. The results, shown in Fig. 6, indicate a lack of reentrant superconductivity effects, at least 
wit hin the range of variation of the control parameters studied here. Special attention was devoted to the region 
{y/U/t, n/U} « {2.5,0.8}, within which low-temperature disorder has been predicted. 19 The numerical simulations 
clearly showed that there were no such effects within this region. In addition, our calculation of the boundary of the 
ordered state using mean-field theory (see Sec. 2) also disagrees with these predictions. Note that Eq. (2), which 
specifies the boundary of the ordered state, is more accurate in this method than is the equation used in the paper 
just cited inasmuch as its derivation did not rely on the assumption that n ;§> 1 [the latter is valid for \ijU ^> 1 (Ref. 
19)]- 

To conclude, we have analyzed the effect of quantum fluctuations in the phase and modulus of the superconducting 
or superfluid order parameter on the character of the ordering in two-dimensional mesoscopic Josephson and granular 
systems within a lattice boson Hubbard model. Quantum Monte Carlo calculations have been used to show that the 
way the system properties change as a result of modulation of the average occupation n umber o f the array elements 
by the chemical potential (the substrate potential) is determined by the parameter q = y/U jtri^ (i.e., the ratio of the 
characteristic Coulomb energy of a granule to the Josephson tunneling energy). For q < 1.5, which is the quasiclassical 
region for the quantum XY model and the region of strong fluctuations in the modulus of the order parameter for 
the Hubbard model (1), the system properties are insensitive to the average number of particles in the granules. In 
the region where there are significant quantum fluctuations in the order parameter (q > 2, 9 < 0.8), we have found 
that the state of the system depends (more distinctly at lower temperatures) on the average number of particles in it. 

This work was supported by grants from the Russian Foundation for Basic Research and the Program on the Physics 
of Solid-Statc Nanostructures. 
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1 - > Two-dimensional Josephson arrays (of mesoscopic elements) with superfluid helium can, in principle, be formed by creating 

the corresponding "picture" in cesium on the substrate (since cesium is not wetted by helium). 1 
2 ) The interference of two Bose condensates has recently been studied. 11 
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